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Abstract We present a hadron production study in the for¬ 
ward rapidity region in ultra-peripheral proton-lead (p + Pb) 
collisions at the LHC and proton-gold (p + Au) collisions 
at RHIC. The present paper is based on the Monte Carlo 
simulations of the interactions of a virtual photon emitted by 
a fast moving nucleus with a proton beam. The simulation 
consists of two stages: the STARLIGHT event generator sim¬ 
ulates the virtual photon flux, which is then coupled to the 
SOPHIA, DPMJET, and PYTHIA event generators for the simu¬ 
lation of particle production. According to these Monte Carlo 
simulations, we find large cross sections for ultra-peripheral 
collisions particle production, especially in the very forward 
region. We show the rapidity distributions for charged and 
neutral particles, and the momentum distributions for neutral 
pions and neutrons at high rapidities. These processes lead 
to substantial background contributions to the investigations 
of collective nuclear effects and spin physics. Finally we pro¬ 
pose a general method to distinguish between proton-nucleus 
{p+A) inelastic interactions and ultra-peripheral collisions 
which implements selection cuts based on charged-particles 
multiplicity at mid-rapidity and/or neutron activity at nega¬ 
tive forward rapidity. 

1 Introduction 

High-energy p+A collisions can be classified into the 
following two categories depending on the impact parameter 
b. In the first category, p+A collisions occur with geomet¬ 
rical overlap of the colliding proton and nucleus, where the 
impact parameter is smaller than the sum of the radii of each 
particle, namely, b <R p +Ra (R p and R,\ are the radius of 
the proton and nucleus, respectively.) 

In the second category instead, the impact parameter 
exceeds the sum of the two radii, b > R p +Ra, thus there is 
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no geometrical overlap between the colliding hadrons and 
hadronic interactions are strongly suppressed. Nevertheless, 
virtual photons emitted from one of the two colliding hadrons 
may anyway interact with another hadron. This process is 
usually referred to as ultra-peripheral collision (UPC, see 
Ref. [1,2] for a review). 

UPCs, so far, have been used for the determination of the 
gluon distribution in protons and nuclei. For example, photo¬ 
production of quarkonium in ultra-peripheral p + A collisions 
can probe a high, or possibly saturated, parton density in pro¬ 
tons at small Bjorken-x (i.e., small parton momentum fraction 
of the momentum of protons). Indeed measurements already 
exist of exclusive J/y/ photoproduction at the CERN Large 
Hadron Collider (LHC), namely, p + Pb —> p + Pb + J / yi [3]. 
Conversely, less attention has been paid, in UPCs, to par¬ 
ticle production in general photon-proton interactions, i.e., 
Y+p —► X, but nevertheless such particle production should 
be considered as well in the investigation of collective nu¬ 
clear effects. Because a large cross section is expected, this 
process in UPCs provides significant background events to 
pure p +A inelastic interaction events (hereafter “hadronic 
interaction”, unless otherwise noted) used for such investiga¬ 
tions. Indeed, a sizable cross section was found for hadron 
production in ultra-peripheral d + Au collisions [4], which 
amounted to ~ 10% of the d + Au inelastic cross section. 
However, in Ref. [4], only the cross section for UPCs was 
presented, and the discussion of the rapidity and momentum 
distributions of the UPC induced events was unfortunately 
neglected. 

In this paper, we discuss the effects of particle production 
by y+p interaction in ultra-peripheral p+A collisions com¬ 
pared to the measurements of hadronic interactions in terms 
of the rapidity and momentum distributions, especially in 
forward rapidity regions at the LHC and the BNL Relativistic 
Heavy Ion Collider (RHIC). Concerning p + Pb collisions 
at y/sm = 5.02TeV at the LHC, we perform the calcula- 
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tions assuming that the measurements of /r°s and neutrons 
are made with zero-degree calorimeters (ZDCs, for example, 
the ATLAS-ZDCs [5]) and the LHCf detector [6], which 
are capable of investigating nuclear effects using hadronic 
interaction events in the very forward region. For the case 
at RHIC, we consider the n° and neutron measurements in 
p + Au collisions at v /siviv = 200GeV in the year 2015. The 
STAR and PHENIX experiments propose a study on partonic 
processes in nuclei using forward prompt photons, where 
decay photons from n°s (from both the hadronic interaction 
and UPCs) would be the dominant background events [7, 
8], Furthermore, measurements of the hadronic-interaction- 
induced prompt photons and 7T°s in transversely polarized 
p + Au collisions may provide key information on the yet 
unestablished contributions of Sivers and Collins effects to 
the single spin asymmetry [7,8], 

Our quantitative discussions on forward hadron produc¬ 
tion are based on Monte Carlo (MC) simulations. The MC 
simulation for UPCs consists of two steps; the virtual photon 
flux is simulated by the STARLIGHT event generator [9, 10] 
and then the subsequent particle production in y+p interac¬ 
tions is simulated by the SOPHIA [11,12], DPMJET [13,14], 
and PYTHIA [15,16] event generators. STARLIGHT in this 
study has been partially customized in order to transfer the 
information on the simulated virtual photon to SOPHIA. The 
MC simulation for hadronic interactions is performed by the 
DPMJET alone. 

The paper is organized as follows. First, in Sect. 2, the 
methodology of the MC simulations is explained. Next, in 
Sect. 3, we discuss the simulation results in terms of the 
rapidity and momentum distributions, where the hadron pro¬ 
duction in UPCs is compared to that in hadronic interactions. 
Additionally, we attempt a reduction in UPC events by re¬ 
quiring associated particles. Conclusions are drawn in the 
last section. In this paper natural units h = c = 1 are used 
throughout. 


2 Monte Carlo simulations Methodology 

As stated above, the MC simulation for UPCs in this 
study consists of two steps. First, we simulate the virtual 
photon flux as a function of the photon energy and impact 
parameter by using STARLIGHT [9], Next, the simulation of 
the y+p interaction is performed by using SOPHIA [11] at 
low energy and either DPMJET [13] or PYTHIA [15] at high 
energy. The methodology of the MC simulation for UPCs 
is explained in the following subsections from 2.1 to 2.3. 
In these subsections, the proton rest frame is referenced to 
unless otherwise noted. The MC simulation for hadronic 
interactions is simply performed by using DPMJET, and will 
be described in Sect. 2.4. 


2.1 Virtual photon flux simulation 


In this paper, the energy spectrum of the virtual photons 
emitted by the relativistic nucleus follows the Weizsacker- 
Williams approximation [17,18] implemented in STARLIGHT. 
The double differential photon flux due to the fast moving 
nucleus with velocity j3 is written as 


d 3 N 
dEydb 2 


Z 2 a x 2 
n 2 Eyb 2 




(1) 


where N is the number of the emitted photons, Ey is the 
photon energy, Z is the electric charge (Z = 82 for Pb and 
Z = 79 for Au), a is the fine structure constant, x = Eyb/y 

__ 1 j2 

(7= \/l-(3 2 is the Lorentz factor), and Kq and K\ are 
the modified Bessel functions. In the case of a relativistic 
nucleus (y^> 1), the contribution of the term K^ix)/y 2 in 
Eq. (1) can be safely disregarded, and in fact STARLIGHT 
considers only the term K 2 (x). For heavy nuclei with a large 
radius, the virtuality of the photon \q 2 \ < (1 /Ra) 2 can be 
neglected. Thus all photons are treated as real photons in the 
simulation for this analysis. Another approximation is due to 
the fact that here we assume a point charge for the nucleus 
and this assumptions may lead to a certain level of systematic 
uncertainty. For example, as discussed in Ref. [19], the pho¬ 
ton flux in reality depends on the choice of the form factor in 
the nucleus by < 20%. 

The probability FuPC(y+p ,x)(b) f° r a single photon in¬ 
teraction with a proton in UPCs as a function of b is given 
by 


[ E T X d 3 N _ 

Ajpc (r+p^x) 0 ) = J Emin dE y d b i °r+p^x {E y )E \ iad (b) dE r 

( 2 ) 

where <3y+ p ^,x ( Ey ) is the total cross section for a single real 
photon interaction with a rest proton and Phad(^) is the prob¬ 
ability of having no hadronic interactions in p+A collisions. 
E'™ n and £™ ax are the minimum and maximum photon ener¬ 
gies. 

In this study, we take Oy+ p ^x{Ey) from the compilation 
of present experimental results [20] when a photon-proton 
center-of-mass energy Wy +p is smaller than 7GeV. A linear 
interpolation is performed between each data point. The cross 
section at the exact photopion production threshold. Ey = 
0.15GeV, for which no experimental measurement exists, is 
forced to zero. At Wy +p larger than 7GeV, (5y+ p -x(Ey) is 
derived from the best COMPETE fit results [20], 

A finite probability for having no hadronic interactions 
P\ rM \(h) is introduced in order to implement a smooth cut off 
for values of the impact parameter approaching: b = R p + Rx- 
Phad(b) is calculated from the Woods-Saxon nuclear density 
and the Glauber model [9], 
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The range of the impact parameter b considered in the 
simulation extends from b mm = 4fm to b max = 10 5 fm. /? mm is 
well below the sum of the effective radii of colliding particles 
(~ 8 fm for both p + Pb and p + Au collisions), and l] VM [(b) 
rapidly approaches zero below 8 fm. The photon energy Ey 
in the simulation ranges from slightly above the photopion 
production threshold, i.e., £™ n = 0.16GeV, to £™ ax . £™ ax 
is obtained from y/b mm and amounts to 700TeV for p + Pb 
collisions at LHC and 1.1 TeV for p + Au collisions at RHIC. 

2.2 Simulation of the low-energy photon-proton interaction 

The particle production from the interaction of a low- 
energy photon with a proton is simulated by the SOPHIA 
2.1 event generator [11], In SOPHIA, particle production via 
baryon resonances, direct pion production, and multiparticle 
production are taken into account. For the baryon resonances, 
the known resonances from A (1232) to A (1950) are consid¬ 
ered with their physical parameters. The resonance decays 
isotropically, depending on the available phase space. The 
non-diffractive interaction, which is implemented using the 
dual parton model [21], starts to dominate at Ey > 2GeV 
increasing with energy. The diffractive interaction is imple¬ 
mented as the quasi-elastic exchange of a reggeon or pomeron 
between virtual hadronic states of the photon and the pro¬ 
ton. The SOPHIA generator is used for the UPC simulations, 
with the photon energy ranging from E™ m to Ey Ut . Ey M is a 
"technical" cut off energy that distinguishes low energy from 
high-energy interactions for SOPHIA and the other generators. 

Here we emphasize that, firstly, the simulation with the 
photon energy Ey < 0.5 GeV is crucial for producing low 
transverse-momentum (p-y) UPC induced events that are 
dominant in the very forward regions of the detector ref¬ 
erence frame (explained in Sect. 3.3), secondly, SOPHIA can 
simulate the interaction of such a low-energy photon with a 
proton above the photopion production threshold, and, finally, 
a newly developed interface to SOPHIA has been introduced 
into STARLIGHT that was not originally coupled to SOPHIA. 
This interface provides two main functions: first, the informa¬ 
tion on the simulated photon by STARLIGHT is transferred to 
SOPHIA, and second, the information on the produced parti¬ 
cles after the simulation of the y+ p interaction are returned 
from SOPHIA to STARLIGHT for, e.g., Lorentz boost, listing 
of the produced particles, etc. 


2.3 Simulation of high-energy photon-proton interaction 

At the photon energy Ey > £.) ut , we perform the simula¬ 
tion of y-l- p interactions by using either PYTHIA 6.428 or 
DPMJET 3.05. STARLIGHT has its own interface to both event 
generators. 


Table 1 Summary of the event generators for y + p interactions and 
their cut off energies. “Low-energy” and “High-energy” in the table 
indicate the energy regions IT ™ 111 < Ey < Ey Ut and E) ut < Ey < E™ ax , 
respectively. 




y+p interactions 



Low-energy 

High-energy 


S OPHIA+PYTHIA 

SOPHIA 2.1 

PYTHIA 6.428 

55GeV 

SOPHIA+DPMJET 

SOPHIA 2.1 

DPMJET 3.05 

6GeV 



In PYTHIA [15], the high-energy photon interactions with 
a proton are classified into three different schemes [22]. Di¬ 
rect events describe the bare photon interaction with a parton 
from the proton, typically leading to high pj jets. In vector 
meson dominance (VMD) events, the photon fluctuates into 
a vector meson and then the vector meson interacts with the 
proton. This class includes low-pT events. Finally, general¬ 
ized VMD events are where the photon fluctuates into a qq 
pair which interacts with a parton from the proton. Single pho¬ 
ton dissociation and single proton dissociation occur in the 
relatively low p-y region. PYTHIA requires a simulated event 
that has a center-of-mass energy Wy +p larger than lOGeV. 
This energy corresponds to a photon energy Ey = 55GeV. 
Thus SOPHIA and PYTHIA are employed for the simulation 
of a y+p interaction for photon energies below and above 
^.cut ,pythia _ Qe V, respectively. 

DPMJET [13] is based on the two-component dual par- 
ton model, y + p interactions are especially implemented 
in the PHOJET MC event generator [23] inside DPMJET. In 
PHOJET, the physical photon is described as a superposition 
of the bare photon and virtual hadronic photon. The bare 
photon directly interacts with partons from the proton. The 
virtual hadronic photon first fluctuates into a qq pair and then 
hadronically interacts with the proton. Both single photon 
dissociation and single proton dissociation are also taken into 
account. DPMJET requires £ , ™ t DPMJET > 6GeV, the lowest 
energy that guarantees usable results from the model. Thus, 
SOPHIA and DPMJET are employed for the simulation of a 
7 +p interaction for photon energies below and above 6GeV, 
respectively. 

As summarized in Table 1, we have thus two types of 
UPC simulations: the first one given by the simulations of 
photohadron production performed by SOPHIA and PYTHIA 
with a cut off energy of /A ut PYTHIA — 55GeV, the second 
one deriving from simulations performed by SOPHIA and 

_ t~' CUt,DPMJET z' i r 

DPMJET with Ey = 6GeV. 

2.4 Simulation of hadronic interactions 

In this paper, DPMJET is used as an event generator for 
the MC simulation of hadronic interactions, which include 
non-diffractive and diffractive interactions but do not include 
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elastic scattering. The multiple scattering process in the inter¬ 
action with a nuclear target, which causes nuclear shadowing, 
is described by the Gribov-Glauber model [24,25] in terms 
of the multiple pomeron exchange. Some of the parameters 
for soft particle production are set at the values that best re¬ 
produce experimental results. The integrated interface CRMC 
1.5.3 [26] is used to access the DPMJET generator. 


3 Predictions for ultra-peripheral collisions at LHC and 
RHIC 

3.1 Total cross sections 


The total cross section for UPCs (<7uPC(y+/>^A')) ca i cu_ 
lated by integrating Eq. (2) over the parameter b: 


rb max 

ff UP C(y+p-+X) = J h<nm P\)?C (y+p^X)db~ 

rb m “ 

= 2k Oy+p-^xiEj) 

Jb mm Je™"' dEydb- 

x Phad(b)bdbdEy. 


(3) 


As discussed in Sect. 2.1, Eq. (1) and l\ a d(b) are ob¬ 
tained by using STARLIGHT, while Gy\ f , ,x(Ey) is taken 
from experimental measurements and best fit results to these 
measurements; thus, Eq. (3) is independent from the other 
event generators: SOPHIA, DPMJET, and PYTHIA. The calcu¬ 
lated cross sections Cupcfy+p^jr) at the LHC and RHIC are 
summarized in Table 2. The cross sections for p+A inelastic 
interactions, calculated by using DPMJET, are also presented. 
We find a sizable ( 7 upc(y+p- 4 X) s > which amount to 20% and 
9% of the hadronic cross sections at the LHC and RHIC 
respectively. Effective cross sections that require forward 71° 
and neutron tagging will be discussed later in Sect. 3.3. 

For reference we have computed the cross section in 
ultra-peripheral d + Au collisions at yts/w = 200 GeV in 
the same manner as Eq. (3) and using the same photon en¬ 
ergy and impact parameter ranges as in Ref. [4], obtaining a 
°UPC(y-H/->x) = 230 mb, which is compatible with the result 
obtained by Ref. [4], Since the photon energy and impact 
parameter ranges used for the comparison are narrower than 
the ranges used in the computation of Table 2, the cross 
section <7[jpc:y ■ d-,X) is smaller than twice the cross section 
^upcfy+p-ix) in Table 2, namely, 2a UPC(y+p ^ x) = 340mb. 
The factor of 2 comes from Oy +( t^>x{Ey) ~ 2c Jy +p ^x(Ey). 

Now, to have a comparison with <7upc(y+p^x)> we have 
also calculated the cross section OuPC(y+A->x) f° r UPCs 
where a photon emitted by a fast moving proton interacts 
with a nucleus A (A is either Pb or Au), although this process 
is not taken into account in the MC simulations performed 
for this study. In ultra-peripheral p +A collisions, the number 
of photons emitted by the proton is generally smaller than 


that by the nucleus due to the Z 2 dependence. The interac¬ 
tions of a single real photon with a nucleus can be roughly 
classified into two categories in terms of the photon energy in 
the nucleus rest frame, namely below or above the photopion 
production threshold. For 7 Me V < Ey < 140 Me V, thus be¬ 
low the photopion production threshold, the cross sections 
for photonuclear absorption processes are studied for Pb and 
Au nuclei in Ref. [27,28]. Here a certain number of neutrons 
are emitted from the decay of the photoexcited nucleus in 
the event, in particular at least one neutron is emitted almost 
100% of the time. For Ey > 140MeV, the y+A cross sec¬ 
tion can be calculated within the Gribov-Glauber approxima¬ 
tion [29]. The calculated y+A cross section Oy+A->x{Ey) is 
compatible with or smaller than the simply scaled y+p cross 
section (Z + N)Oy+ p ^,x{Ey) (Z and N are the atomic num¬ 
ber and neutron number of the nucleus, respectively.) The 
suppression of (£ 7 ) relative to (Z + N)<Jy+ p ^,x(Ey) 

appears at Ey > 2 GeV because of nuclear shadowing [29]. 
Consequently, in ultra-peripheral p + Au and p + Pb colli¬ 
sions, the process with photon emission from the proton is 
suppressed with respect to that from the nucleus. The cross 
section of the former process (C7upc(y+A->X)) amounts to 6 % 
of the latter (CT UPC{r+p ^ x) ). 

Since the UPC simulations used in this paper assume that 
only a single photon is produced from the moving nucleus 
in any event, UPCs involving two or more photons are not 
taken into account in the study. However UPC processes with 
the exchange of two photons between the proton and the 
nucleus generally have a huge cross section. In particular, 
the two-photon exchange leading to di-electrons has a cross 
section 0 \;i>c!yi y- ,,e-e) °f ~ 29b at the LHC and of ~ 4b 
at RHIC (oc (Z p a) 2 (Z A ct) 2 , where Z p and Z,\ are the elec¬ 
tric charges of the proton and nucleus, respectively). These 
cross sections are obtained by scaling the corresponding 
cross sections in ultra-peripheral Pb+Pb and Au+Au colli¬ 
sions [30,31] with the ratio of (Z p /Za) 2 , respectively. The 
two-photon exchange leading to other particle pairs, e.g., 
p + p~, T + T~, and mesons, are at most 10 ^ 3 % compared 
with the di-electron channel [32]. Concerning forward 71° 
and neutron productions providing background contribution 
to hadronic interaction events, we consider a higher-order 
process where the two-photon exchange is accompanied by 
an additional inelastic interaction of the single photon emit¬ 
ted from the nucleus with the proton. This process can be 
factorized into two subprocesses: the two-photon exchange 
and the y+p interaction [33]. Accordingly, the cross section 
is calculated using: 

®UPC(y+y->e+e, y+p->X) = 

fb max 2 (4) 

J b mm ^ UPC (r+r^e+e) (b)Pupc(y+p—>x) ( b)db 

where P\jpc(y+y^>e+e) (P) is the probability for the two-photon 
exchange leading to di-electron and PuPC(y+p-D<:) P) is the 
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Table 2 Cross sections for particle production in ultra-peripheral colli¬ 
sions and hadronic interactions at the LHC and RHIC. 




UPC (mb) 

Hadronic interaction (mb) 


Total 

7T 

n 

Inelastic 


n 

LHC 

434 

78 

153 

2189 

91 

125 

RHIC 

170 

9 

73 

1851 

67 

35 



probability defined in Eq. (2). With the former probabil¬ 
ity taken from the results in Refs. [34,35], we obtain 
®UPC(y \ y--,e-e. yip -X) = 66 mb at the LHC and 15 mb at 
RHIC. These cross sections correspond to 15 % and 9% com¬ 
pared with OuPC(y+p->x) shown in Table 2 at the LHC and 
RHIC, respectively. 

3.2 Rapidity distributions 

The charged and neutral particle pseudorapidity (t7i a b) 
distributions in the detector reference frame are shown in 
Fig. 1. The solid curves indicate the UPC simulation events 
generated by using SOPHIA at low energy and DPMJET at 
high energy. The dashed curves indicate the UPC simulation 
events generated by using SOPHIA at low energy and PYTHIA 
at high energy. The dotted curves indicate the predictions 
for p + Pb and p + Au inelastic events with DPMJET at the 
LHC and RHIC respectively. Hereafter the directions for the 
moving proton and nucleus are assumed to have positive 
and negative rapidities respectively. The directions for both 
particles are indicated by the arrows in the upper left panel 
of Fig. 1 (p for proton and A for nucleus). 

In Fig. 1, the pseudorapidity distributions in hadronic 
interactions achieve a plateau at mid-rapidity and also have a 
large number of spectator nucleons at rji a b ~ —8. Conversely, 
particle production in UPCs is clearly clustered around posi¬ 
tive forward rapidities, since the UPC induced events in this 
study are produced only by the interactions of the photon 
emitted from the nucleus with the proton. It should be noted 
that the cross sections for UPCs exceed those of hadronic 
inelastic interactions for charged and neutral particles at 
t]i a b > 9 at the LHC and rji a b > 7 at RHIC. This indicates that 
contamination from background UPC events, could spoil the 
investigation of collective nuclear effects and spin physics 
carried out from measurements of hadronic interactions in 
those rapidity regions. 

3.3 Transverse and longitudinal momentum distributions 

The simulated p-y and longitudinal momentum fraction 
( z , defined as p v j p, max ) distributions for 7T°s and neutrons 
at positive forward rapidities (direction of the proton rem¬ 
nant), are shown in Fig. 2 for the LHC and in Fig. 3 for 
RHIC. For the distributions at the LHC, we chose rapidity 



-10 0 10 -10 0 10 
Pseudorapidity Pseudorapidity 


Fig. 1 Charged (left) and neutral particle (right) pseudorapidity distribu¬ 
tions at the LHC (top) and RHIC (bottom), respectively. The solid curves 
and dashed curves indicate the UPC simulation events generated by us¬ 
ing STARLIGHT+SOPHIA+DPMJET and STARLIGHT+SOPHIA+PYTHIA, 
respectively. The dotted curves indicate the simulated p + Pb and p + Au 
inelastic events with DPMJET at the LHC and RHIC, respectively. The 
directions of the moving proton (p) and the nucleus (A) are indicated 
by arrows in the upper left panel. 

regions in the detector reference frame 8.5 < ji ab < 11.0 for 
K °s and 7.0 < yi a b < 9.05 for neutrons, which roughly corre¬ 
spond to the acceptances of the ATLAS-ZDCs [5] and LHCf 
detector [6]. These detectors provide an opportunity to inves¬ 
tigate the effects of high parton density on forward 7T°s and 
neutrons, which emerge as a suppression in the momentum 
distributions in p + Pb inelastic interactions relative to that 
of p + p inelastic interactions. 

In the upper left and bottom left panels of Fig. 2, the 
P'\ distributions of /r°s and neutrons in UPCs show a steep 
peak at py S3 0.2 GeV. These peaks originate from the 
Y+p —► 7T° +p and Y + P —> K + + n channels, via baryon 
resonances. In fact, in the proton rest frame the y+ p inter¬ 
actions, with a photon energy ranging from F™ 111 to 0.5 GeV, 
have a center-of-mass energy of 1.1 < Wy +p < 1.3 GeV and 
thus occur in the baryon resonance region, which has a larger 
cross section compared to other energy regions. Conversely, 
the y + p interactions with higher photon energies are sup¬ 
pressed due to a decrease in the photon flux with increasing 
photon energy. Therefore the n °s and neutrons emitted by 
the decay of the baryon resonances due to low-energy Y+P 
interactions (dominantly A + (1232)) which typically have 
p-y S3 0.2 GeV provide substantial contributions to the p-y 
distributions of each particle. Thus, the double differential 
UPC cross sections exceed those of hadronic interactions 
for py- ss 0.2GeV. The dominance of the y + p —» n + + n 
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channel in UPCs is also evident in the bottom right panel of 
Fig. 2. Forward neutrons produced in UPCs have a larger z 
value; low momentum neutrons produced by a low-energy 
7 +p interaction in the proton rest frame are boosted to nearly 
the same velocity of the projectile proton. Finally, we see 
that the presence of UPCs certainly provides a significant 
background contribution to the study of collective nuclear 
effects. 

The pj and z distributions for 7T°s and neutrons produced 
at RHIC are shown in Fig. 3. We chose the rapidity regions 
3-1 < yiab < 3.8 for 7T°s and 4.0 < yi a b < 5.4 for neutrons, 
which correspond to the acceptance of the MPC-EX detec¬ 
tor [8] and ZDC [36] of the PHENIX experiment respectively. 
Since these rapidity regions are near the acceptance of the 
FMS+pre-shower detector [7] and ZDC [36] of the STAR 
experiment, the following discussion is essentially applicable 
also to measurements performed at the STAR experiment. 
The prompt photon measurements in p + Au collisions with 
the MPC-EX detector in PHENIX and the FMS+pre-shower 
detector in STAR could provide key information on par- 
tonic processes in the Au nucleus, whereas photons from 
7T°s decays, produced both from hadronic interactions and 
UPCs, could be the dominant background events. The pj 
distributions of 7T°s from UPCs have almost the same shapes 
as those coming from hadronic interactions (the upper left 
panel), while the absolute yield is at most 3% of that of 
hadronic interactions. Thus, we conclude that UPCs provide 
a negligible contribution to the amount of n°s giving back¬ 
ground photon events. The z distribution of 7T°s from UPCs 
(the upper right panel) is slightly steeper than for hadronic 
interactions, but the absolute yield is once again negligible. 
On the other hand, neutrons produced in UPCs compete with 
those from hadronic interactions at pj 0.2GeV. This can 
be explained by the same mechanism found in Fig. 2, namely 
dominance of the y +p -+ n + + n channel via baryon reso¬ 
nances in UPCs. The z distribution of neutrons has a similar 
shape to that shown in the bottom right panel of Fig. 2. Con¬ 
cerning the measurements of forward prompt photons, 7t° s, 
and neutrons in polarized p + Au collisions, which are sen¬ 
sitive to the origin of spin asymmetry [37], we confirm that 
the UPC contribution to the total number of prompt photons 
and 7T°s is negligible, while neutrons instead, are produced in 
similar amounts from both UPCs and hadronic interactions. 

We need to test now, both for the LHC and RHIC simu¬ 
lations, a possible dependence of the pj and z distributions 
on the cut off energy Ey Ut that was introduced in Sect. 2.2 
and 2.3. The comparison of the two distributions at the LHC, 
with one based on the UPC simulations with DPMJET above 
EyUt DPMJET _ 6Q e y anc j t | lc other based on the UPC simu¬ 
lations with DPMJET above £™ tJJPV1JEI = 55 GeV, shows a 
negligible difference between these two. The comparison in 
the RHIC case shows larger differences than at LHC; never¬ 
theless they are at the same level of those between DPMJET 


and PYTHIA. There is also a significant difference for the pj 
distributions of neutrons at the LHC between DPMJET and 
PYTHIA (the bottom left panel of Fig. 2); PYTHIA predicts a 
harder spectrum and has ~ 8 times larger number of neutrons 
at pj ~ 2 GeV than DPMJET. This difference is mostly caused 
by a strong dependence of the multiplicity of leading baryons 
on the minimum pj (Pxmin) for the multiple interactions (i.e., 
hard scattering) implemented by each model. In general, a 
lower pj min provides more multiple interactions in an event 
which then leads to a higher multiplicity. Changing the pxmin 
value in PYTHIA (PARP(81), 1.9GeV as a default) to the 
default value of DPMJET (2.5GeV) significantly modifies 
the pj distribution in PYTHIA bringing an overall agreement 
between the two. 

Looking back at Table 2 we see that the effective cross 
sections for n {) and neutron productions are defined as the 
cross sections having at least one K° or neutron hitting within 
the rapidity ranges highlighted above. We find that UPC cross 
sections are similar or larger than the hadronic cross sections, 
except for the n° one at RHIC. It should be noted that the 
effective cross section estimation involves simulations of 
specific hadron productions and is thus no longer independent 
of the generators used for hadronic interactions. The values 
in Table 2 for UPCs are calculated by using STARLIGHT, 
SOPHIA, and DPMJET, and those for hadronic interactions 
are calculated by using DPMJET. Nevertheless, the essence 
of our conclusions will not depend on the choice of event 
generators. 

Finally, we roughly estimate the contribution to UPCs 
where the proton acts as the photon source, which was not 
taken into account in the MC simulation of this paper, and the 
possible changes to the pj and z distributions. Because of the 
directions of the colliding photon and nucleus and because 
of the smaller cross sections (calculated in Sect. 3.1), the 
number of 7T°s and neutrons produced at yi a b > 3 is negligible 
compared with that from UPCs where the nucleus acts as the 
photon source. Conversely, the two-photon exchange process 
followed by the y+p interaction, currently not implemented 
in the MC simulation as well, would yield n°s and neutrons 
at positive forward rapidities amounting to 15% of those 
shown in Fig. 2 at the LHC and to 9 % of those shown in 
Fig. 3 at RHIC. A more detailed study would need to take 
two-photon and di-electron interactions into account in the 
MC simulation framework and is thus beyond the scope of 
the present paper. 

3.4 Reduction of ultra-peripheral collisions contributions 

As seen in Fig. 1, particle production in UPCs is clus¬ 
tered at positive forward rapidity regions, whereas those in 
hadronic interactions show a plateau at mid-rapidity and have 
also a large number of spectator nucleons at T)i a b ~ —8. In 
this section, we present the two methods to separate UPCs 
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Fig. 2 Simulated px and z spectra for 7t°s and neutrons in p + Pb collisions at the LHC. The solid curves and dashed curves indicate the UPC 
simulation events generated by using STARLIGHT+SOPHIA+DPMJET and STARLIGHT+SOPHIA+PYTHIA, respectively. The dotted curves indicate 
the simulated p + Pb inelastic events with DPMJET. 





Fig. 3 Simulated px and z spectra for TT s and neutrons in p +Au collisions at RHIC. The solid curves and dashed curves indicate the UPC 
simulation events generated by using STARLIGHT+SOPHIA+DPMJET and STARLIGHT+SOPHIA+PYTHIA, respectively. The dotted curves indicate 
the simulated p + Au inelastic events with DPMJET. 
















































and hadronic interactions by exploiting those differences in 
the rapidity distributions. 

First, we investigate the effects of requiring some activi¬ 
ties in the mid-rapidity region. The following cuts are applied 
to inclusive n° and neutron measurements to eliminate UPC 
induced events as well as to ensure that the hadronic inter¬ 
action events remain unchanged: (1) the number of charged 
particles should be greater than 2, (2) the charged particles 
should have pj > 0.2GeV, and (3) the charged particles 
should have rapidity |rji a b| < 2.5 at the LHC and | Tji a b | < 0.35 
at RHIC. The rapidity regions used in the cuts correspond 
to the rapidity ranges of the ATLAS and PHENIX inner 
detectors [38,39]. 

Figure 4 shows the py distributions after the cuts at the 
LHC and RHIC. The absolute yields of UPCs at the LHC 
(RHIC) are reduced to less than 60% (8%) for 7T°s and 
40% (10%) for neutrons, even though the absolute yields 
of hadronic interactions are kept to be larger than, at most, 
85% (50%) for 7T°’s and 85% (20%) for neutrons. Thus 
these cuts on the charged particles at mid-rapidity reduce 
the relative yields of UPCs to hadronic interactions by, at 
most, 15% (0.2%) for 7T°s and 25% (1.5%) for neutrons. 
However the cuts still leave a remnant contribution from 
UPCs and unavoidably reduce hadronic interaction events. 
It should also be noted that the rejected hadronic interaction 
events are mostly single and double diffractive events, gener¬ 
ally characterized by a small number of charged particles at 
mid-rapidity, and that the reduction efficiency of 10% — 50% 
highly depends on forward n {) and neutron energies. 

Next, we test sharper cuts, which require an activity in 
the negative very forward region, i.e., the direction of the 
nucleus remnant, and which can be tagged by a ZDC. In 
this rapidity region, only hadronic interactions produce a 
large number of spectator nucleons fragmented from the 
colliding nucleus. Each nucleon has an approximate energy 
of 1.58TeV at the LHC and lOOGeV at RHIC. The cuts 
applied to the simulated events consist of three requirements: 
(1) the number of neutrons should be greater than 1 (note 
that a proton can be swept away by the magnets located 
between an interaction point and a detection point, and thus 
no proton reaches the detector), (2) the neutrons should have 
E > 1 TeV at the LHC and E > 50 GeV at RHIC, and (3) the 
neutrons should have a rapidity r)] a b < —6.5. The rapidity 
region r)i a b < —6.5 roughly corresponds to the acceptance of 
the ZDC located in the nucleus-going side. As was inferred 
in Fig. 1, we find that the contribution of UPCs is efficiently 
eliminated by these cuts, keeping the number of hadronic 
interactions unchanged. 

We conclude that a reduction in the contributions from 
UPCs to the measurements of forward hadrons is certainly 
feasible by requiring some activity in mid and forward ra¬ 
pidity regions. In particular, we expect a strong reduction if 
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Fig. 4 Simulated px spectra for 7t°s and neutrons with the cuts ap¬ 
plied at the LHC and RHIC (see text in detail). The solid curves and 
dashed curves indicate the UPC simulation events generated by us¬ 
ing STARLIGHT+SOPHIA+DPMJET and STARLIGHT+SOPHIA+PYTHIA, 
respectively. The dotted curves indicate the simulated p+A inelastic 
events with DPMJET. 


we detect spectator nucleons, for example, with a ZDC, at 
negative rapidity rji a b < —6.5. 


4 Conclusions 

Hadron production in forward rapidity regions by ultra¬ 
peripheral p + Pb collisions at the LHC and p + Au colli¬ 
sions at RHIC was discussed. The present paper was based 
on MC simulations of the interaction of virtual photons emit¬ 
ted by a fast moving nucleus with a proton, by using sev¬ 
eral event generators: STARLIGHT, SOPHIA, DPMJET, and 
PYTHIA. STARLIGHT in this paper was customized in order 
to transfer the information on the simulated virtual photon to 
SOPHIA and to receive the information on the produced parti¬ 
cles from SOPHIA. This modification made possible the sim¬ 
ulation of UPC induced events starting from the photopion 
production threshold. We found large cross sections for 7T° 
and neutron productions in the very forward region, leading 
to a substantial background contribution to the measurements 
of hadronic interactions at both the LHC and RHIC. There¬ 
fore, the presence of UPCs had to be taken into account in the 
analyses focused on the investigation of collective nuclear 
effects and spin physics, in order to correctly evaluate the 
fraction of the hadronic cross section relative to the mea¬ 
sured cross section. We propose two types of cuts to reduce 
the fraction of UPCs relative to hadronic interactions; one 
requires the presence of charged particles at mid-rapidity. 
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while the other requires the presence of spectator neutrons 
at negative very forward rapidity. The former cut certainly 
reduces the UPC induced events while unavoidably rejecting 
some hadronic events. The second cut efficiently eliminates a 
fraction of forward n°s and neutrons produced by UPCs and 
does not change the number of hadronic interaction events. 
The proposed methods in this paper were simply based on 
a tracking detector and a ZDC and thus are generally appli¬ 
cable to the measurements in other experiments that have a 
similar detector design. 
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